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We study the physical properties of a ballistic heterostructure made of a ferromagnet (FM) and 
a spin-triplet superconductor (TSC) with a layered structure stacking along the direction perpen- 
dicular to the planes where a chiral p x +ip y pairing occurs and assuming spin dependent processes 
at the interface. We use a self-consistent Bogoliubov-de Gennes approach on a three-dimensional 
lattice to obtain the spatial profiles of the pairing amplitude and the magnetization. We find that, 
depending on the strength of the ferromagnetic exchange field, the ground state of the system can 
have two distinct configurations with a parallel or anti-parallel collinearity between the magnetic 
moments in the bulk and at the interface. A magnetic state having non coplanar interface, bulk, 
and Cooper pairs spins may be stabilized if the bulk magnetization is assumed to be fixed along a 
given direction. The study of the density of states reveals that the electronic reconstruction in the 
FM plays an important role in the setting of the optimal magnetic configuration. We finally discuss 
the existence and the character of induced equal-spin pair correlations in the FM-TSC system. 

PACS numbers: 74.45.+C, 74.20.Rp, 74.50.+r 



I. INTRODUCTION 

The past few decades have been marked by a grow- 
ing interest in the study of the interplay between 
superconductivity and ferromagnetism in heterostruc- 
tures both for the potential application in the field of 
spintronics and because of the underlying fundamen- 
tal physics^— . The control of the spin degree of free- 
dom in charge currents that can flow without dissipa- 
tion represents one of the major challenges in the con- 
text of hybrids based on superconductors and ferromag- 
nets. In this respect, the presence of spin-triplet pair cor- 
relations in ferromagnet(FM)-superconductor(SC) het- 
erostructures is expected to play a key role in the genera- 
tion of new types of spin- valves or unconventional Joseph- 
son junctions. 

One of the major obstacle to such achievement is the 
difficulty in obtaining a coexistence of superconductiv- 
ity and ferromagnetism within a uniform quantum state. 
The competition between the two orders can be reduced 
by a phase or amplitude modulation of the superconduc- 
tor order parameter, as for the so called Fuldc-Ferrell- 
Larkin-Ovchinnikov (FFLO) stated, by having ferromag- 
netic domains with no net magnetization in average on 
the scale of the superconducting coherence length, or via 
specific mechanisms driving the ferromagnetism^. Such 
incompatibility, due to the opposite nature of the order 
parameters, gets less severe in a ferromagnet- spin-singlet 
superconductor(SSC) heterostructure through the occur- 
rence of nanoscalc proximity effects at the interfaces. One 
remarkable feature due to the proximity effects is an in- 
duced pairing amplitude in the ferromagnet, character- 
ized by a damped oscillatory behaviour— ~—. This leads 
to a nonmonotonous dependence of the superconducting 
critical temperature as well as a 0-ir transition in the 
Josephson current and the possibility to suppress the su- 



perconductivity according to the FM thickness 2 - ' 10 ^ 1 ' 13 . 
Moreover, spin-triplet pair correlations being even in 
space and odd in time or frequency can also be induced 
in the presence of magnetic inhomogeneities^— . Non 
uniform magnetic configurations can be intrinsic or arti- 
ficially designed, as for instance in the case of a FM1- 
FM2-SSC heterostructure with noncollinar magnetiza- 
tion between FM1 and FM2 ferromagnets. The FM2 
ferromagnet acts in this configuration as a spin-active in- 
terface where its magnetization orientation with respect 
to the FM1 magnetization is the key parameter to tune 
the proximity effects. In turn the presence of the FM2 
layer results into a change of the pair-breaking intensity 
and of the superconducting critical temperature. Both 
the proximity effects and the induced spin-triplet pairing 
can contribute to generate a spin-valve effect by turn- 
ing the superconductivity on and off depending on the 
pattern of the magnetic profile^—. 

Differently from the induced triplet correlations, a 
spin-triplet superconductor (TSC) is marked by a bro- 
ken symmetry configuration due to a characteristic pair- 
ing potential. Spin-triplet superconducting order has 
been predicted in systems like the heavy-fermion com- 
pounds UPt^ and UGe^ 3 -, the quasi-one-dimensional 
organic Bechgaard salts (TMTSF) 2 X^i, and the per- 
ovskite ruthenium oxide Sr2RuO,f2£. While in the first 
two cases a non-unitary pairing seems to occur, a uni- 
tary chiral-state where time-reversal symmetry is bro- 
ken by orbital degrees of freedom is realized in the lat- 
ter—. A peculiarity of this state is its non-trivial topo- 
logical character, specified by a fully open excitation en- 
ergy gap coexisting with topologically protected gapless 
edge states. Spin-triplet pairing may also manifest itself 
in non-centrosymmetric superconductors (NCS) 27 i 28 to- 
gether with the spin-singlet component, due to the lack 
of a crystal inversion center that leads to an asymmet- 
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FIG. 1. (color online) FM-TSC heterostructure with the 
interface perpendicular to the z-direction. The blue layers 
are for the spin-triplet superconductor with a d-vector (thick 
black arrow) perpendicular to the (sy)-plane and the red lay- 
ers indicate the ferromagnetic region. The red thin arrow 
stands for the magnetization M in the interior of the ferro- 
magnetic region (at the interface) forming a relative angle 9 
(p) with respect to the d-vector. L x and L y are the in-plane 
lateral size of the heterostructure. 



ric spin-orbit interaction and in turn results into parity 
mixing^. The increasing number of materials having a 
superconducting order parameter with a spin-triplet sym- 
metry makes more achievable the design of TSC based 
heterostructures. In particular the recent preparation 3 ^ 
of superconducting thin films of S^RuO^ 3 -^ provides a 
key milestone for the development of such hybrid sys- 
tems. Within this context the epitaxial growth of the cu- 
bic itinerant ferromagnct SrRu03 on the top of Si'2Ru04 
could allow one to obtain an FM-TSC heterostructure 
with the desired control of the interface. 

The physical behaviour of the FM-TSC heterostruc- 
ture is expected to be quite different from that of the FM- 
SSC. This is especially due to the basic different structure 
of the spin-triplet order parameter that has a vectorial 
nature and thus can couple to the orientation of the FM 
moment as well as to the orbital one. Such features play 
a distinct role in controlling the dependence of the prox- 
imity effects in TSC-FM bilayers with diffusive FM22 and 
the sign of the current in TSC-FM-TSC Joscphson junc- 
tions^. The interface with the FM also induces specific 
effects in the TSC generating a Joscphson-likc coupling 
between the spin t and 4- condensates of the TS C 34 ' 35 . 
Furthermore, in a FM-TSC heterostructure with a planar 
geometry the ferromagnet magnetization tends to favor a 
specific orientation in a way that depends on the orbital 
symmetry of the order parameter rather than that of the 
Cooper pairs spin 3 ^. The role of spin-misalignment on 



hybrids based on TSC, in analogy with the FM1-FM2- 
SSC systems, can significantly modify the above scenario 
and yield novel physical behaviors. So far it has been 
considered only for the case of a half-metal ferromagnct 
demonstrating that extra midgap states can be induced 
as well as a nonzero charge conductance 3 ^. 

In this paper we aim to study the role of spin- dependent 
interface processes in determining the proximity effects 
and the physical properties of an FM-TSC heterostruc- 
ture. In particular, we consider FM and TSC layers 
stacked along the z-direction perpendicular to the plane 
where a chiralp^+ipy pairing occurs in the TSC and with 
a spin-active interface such as to mimic an FM1-FM2- 
TSC heterostructure (see Fig. [J). We employ a fully self- 
consistent Bogoliubov-de Genncs approach on a lattice to 
determine the proximity behaviour through the spatial 
profile of the spin-triplet pair amplitude and the magne- 
tization. We obtain the optimal orientation of the ferro- 
magnet magnetization with respect to the interface-bulk 
magnetization misalignment in the presence of the spin- 
triplet superconductor. The interface-bulk misalignment 
tunes the coupling between ferromagnctism and super- 
conductivity by leading to a subtle competition between 
the energy contributions associated with the proximity 
effects and with the superconducting order parameter 
and those arising from the electronic reconstruction at 
the interface. We demonstrate that, depending on the 
strength of the ferromagnetic exchange, the ground state 
of the system can have two distinct magnetic configura- 
tions with in-plane collinear parallel or antiparallel mag- 
netic moments between the interface and the interior of 
the ferromagnet. The transition between the two config- 
urations, tuned by the strength of the ferromagnetic ex- 
change in the bulk or at the interface, is strongly related 
to the electronic reconstruction occurring in the ferro- 
magnetic region of the system. On the other hand a non 
collinear magnetic moment at the interface emerges with 
respect to a given magnetization in the bulk that is not 
coplanar to the spin of the Cooper pairs. Finally, we dis- 
cuss how equal-spin triplet pair correlations are induced 
at the FM-TSC interface depending on the magnetiza- 
tion orientation in the ferromagnet with respect to the 
direction of the Cooper pairs spin. We believe that the 
study of such a geometric configuration has a relevant 
impact in view of a possible application to a Sr2Ru04 
heterostructure stacked along the c-axis. 

The paper is organized as follows. In Sect. II we intro- 
duce the model Hamiltonian and the relevant physical 
quantities to be determined. The Sect. Ill is devoted 
to the proximity effects in the FM-TSC system while in 
Sect IV we present the results for the FMi-FM 2 -TSC case 
and the phase diagram corresponding to the energetically 
most favourable magnetic configurations. In Sect. V wc 
analyze the induced equal-spin pair correlations at the 
interface and finally the concluding remarks arc reported 
in Sect. VI. 
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II. THE MODEL 

We consider a FM-TSC heterostructure with a lay- 
ered geometry (see Fig. [1]) described by a single-band 
tight-binding model Hamiltonian with an in-plane at- 
tractive nearest-neighbor interaction such as to yield a 
spin-triplet pairing with chiral symmetry in the super- 
conducting side. In the framework of an effective Stoncr 
model for itinerant magnets, an exchange field h (h- m t) is 
introduced to yield a nonzero spin polarization in the fer- 
romagnetic layers (at the interfacial layer with the TSC), 
respectively. Denoting by z and {x, y} the directions 
perpendicular and parallel to the interface layer, respec- 
tively, we assume that the system is uniform along the x 
and y axes, so that the translational symmetry is broken 
only in the z direction. The total Hamiltonian of the 
system is then defined on a lattice with size L x x Ly x Ijj 
(L x = L y = L) with periodic boundary conditions along 
x and y. The simulation is performed for a system size 
with L = 40 and L z = 40. We indicate each site by a 
vector i = (i,i z ), with i labeling the (xy) in-plane atomic 
positions and i z the different layers along the z-direction. 
The spin-active layer is placed at the site i z = between 
the FM region (i z < 0) and the TSC one (i z > 0). The 
hopping amplitude is given by t x = t y = t for the in-plane 
processes and t z = t/2 uniformly for all the inter- layers 
charge transfer, including that at the FM-TSC interface. 
A different choice of the single particle electronic parame- 
ters and the system size alters only quantitatively and not 
qualitatively the presented results. All the energy scales 
are in unit of t and for clarity, the exchange fields have 
been scaled with respect to the amplitude /ihm = 3.15 1 
corresponding to the ferromagnetic half-metallic state, 
i.e. when there is only one type of carriers at the Fermi 
level and the minority spin electrons have a gap in the 
excitation spectrum. With the above assumptions the 
Hamiltonian can be expressed as 

H = Y ( efc(T ~ M)cL(^) C feafe) + 

- Y V [nit(iz)nji(i z ) + n^(i z )n jt (i z )] 

<ij>,iz>0 

+ Y £*«( C L(*« + *) c far(*») + h-c.) + 
<S=±1 ka 



Y. h ■ s(i z ) - hi n t ■ s(0) 



(1) 



iz<0 



with Ckaiiz) being the annihilation operator of an elec- 
tron with planar momentum k, spin a at the layer i z and 
where (ij) denotes the in-plane nearest-neighbor sites, 
(i is the chemical potential, V is the in-planc pairing 
strength, s(i z ) = J2i, s , s < c\ s {i z )a s ^c ls '{i z ) is the spin 
density of the layer i z , and tka = — 2t[cos(k x ) + cos(k y )] 
is the in-planc tight-binding spectrum. 

For the analysis of the proximity effects and of the 
most energetically favorable magnetic configurations the 



quartic term of the model Hamiltonian in Eq. [T] is de- 
coupled in the Hartree-Fock approximation within the 
pairing channel. We introduce the pairing amplitude for 
an in-plane bond Fij(i z ) = (ci^(i z )cj i(i z )} at a given 
layer i z . The average (^4) indicates the thermal expec- 
tation value of the operator A. Within such decoupling 
scheme, the Hartree-Fock HHamiltonian is 



Hhf= ( Ck ° ~ M)cL(^) c fcafe) + 

k,i z ,(T 

+ V [E l3 (2 z )cl(* z )ci t (i z ) 

<ij>,i z >0 

+ FjiiQc^iQc]^) + h.c. 



5=±1 ka 

Y h ■ s(i z ) - hint ■ s(0). 



(2) 



From the pairing amplitude Fij(i z ) it is possible to 
build up the superconducting order parameter with 
the desired symmetry, to be then determined self- 
consistentl y 38 ! 39 . By a proper selection of the spin-triplet 
channel one can derive the pairing amplitude for each 
layer by summing up the contributions over all the in- 
plane k vectors. 

In general, for a spin-triplet superconductor the order 
parameter can be expressed in a matrix form as 



A(k) 



A n (k) A n (k) 
A it (fc) A u (k) 



id,, 



where the (i-vector components are related to the pair 
correlations of the various spin-triplet configurations. In- 
deed, d x = i(-A tt (fc) + A u (fc)), d y = ±(A n {k) + 
A^(k)), and d z = A^(k) are expressed in terms of the 
equal spin A^f(k) and A^(fc), and the anti-aligned spin 
A-\-i(k) pair potentials. For the model upon examina- 
tion the pairing interaction V is non zero in the chan- 
nel and thus the only non- vanishing order parameter is 
A^(k). This implies that the d- vector is oriented along 
the z-dircction, which is chosen to be perpendicular to 
the layers as indicated in Fig. [TJ 

It is worth pointing out that since we are dealing with a 
hybrid structure where the superconductor is interfaced 
to a ferromagnet, the spin-triplet pair correlations can 
have also contributions in the tt an d 4-4- channels. For 
this reason it is useful to introduce the following pair 
correlators 

F pAiz)= y~T~ y" sm ( k a)(ckt(iz)c-]a(iz)) (3) 

lj X Ly k 

F™{i z )= ^sin(k Q )(c kff (i z )c_ k(T (i z )) (4) 

x y k 

with a — x,y and a =t, i- 
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The analysis is performed by choosing the chemical po- 
tential (p = — 1.8t) and the pairing strength (V = — 2.5i) 
such that the most stable superconducting state has 
an orbital symmetry of p x +ip y type. Concerning the 
ferromagnetic side of the heterostructure the inner layers 
magnetization M depends on the exchange field h, and 
can form an angle 9 with respect to the direction of 
the d- vector (Fig. [IJ. On the other hand, the interface 
exchange h- lnt leads to a magnetic moment forming an 
angle p with respect the d- vector (Fig. [T]). Due to the 
rotational spin symmetry of the superconducting state, 
the x- and y- directions are equivalent. Hence, we can 
consider h to lie in the (a;z)-plane with components 
given by h = ft. (sin(0), 0, cos(0)) and similarly for 
hmt = /iint(sin(p),0, cos(p)). 

Then, after having diagonalized the Hartree-Fock 
Hamiltonian -ffjjF (Eq.©), the Gibbs free energy can 
be determined as 

G= ~ In Tr[c-' 3ifHF ] (5) 

where /3 = -j^~r * s ^ ne mverse temperature, kg being the 
Boltzmann constant. In order to separate the energy con- 
tributions resulting from the proximity effects and from 
the modification of the energy spectrum by those that 
enters in the superconducting condensation energy via 
the change of the order parameter in the TSC side of the 
heterostructure, it is useful to introduce the quantity E op 
expressed as 

E ap = -(-p\ [\F P AQ\ 2 + \F Pv (i*)\ 2 ] ■ (6) 



Such quantity reflects the changes of the superconducting 
order parameter in the FM-TSC heterostructure related 
to the magnetization configuration of the ferromagnetic 
layers. In particular, a minimum in E op implies a maxi- 
mum amplitude of the superconducting order parameter. 

Furthermore, for clarity the results are presented by 
rescaling the Gibbs free energy and the order-parameter 
derived energy to their minimum value. Indeed, we in- 
troduce the quantities AG = G - G min and AE'op = 
E op — -E"p ln where G min and E™ n are the minimum am- 
plitude of the Gibbs energy and of the order parameter 
energy, respectively, in the analyzed range of parame- 
ters. This rescaling allows one to directly extract the 
energy scale that sets the variation of the Gibbs energy 
landscape and to compare different energy profiles as a 
function of the physical parameters involved in the anal- 
ysis. Still, the variation of AG and A_E op are presented 
by scaling their amplitude to the value of E op associated 
with a single superconducting layer in the homogenous 
state. 



III. PROXIMITY EFFECTS AND OPTIMAL 
MAGNETIC CONFIGURATION OF FM-TSC 
SYSTEM 

In this section we present the evolution of the pairing 
amplitude of the FM-TSC system assuming that the 
ferromagnet has an uniform magnetization i.e., 9 = p. 
We consider for the discussion a representative value of 
the ferromagnetic exchange h = h lnt = 0.16. 

As depicted in Fig. [2] the spatial dependence of the 
real part of the pairing amplitude for the p x compo- 
nent has a distinct behavior with respect to the angle 
9 when considering the FM and the TSC side of the het- 
erostructure. The imaginary part of the Pj, component 
exhibits the same behaviour of the p x counterpart due 
to symmetry. In the TSC region the pairing amplitude 
is reduced close to the interface, with a maximum re- 
duction for 9 = (see the inset of Fig. [2]). This trend 
is expected when the pairing amplitude variation is con- 
trolled by the Andreev reflections at the interface in such 
a way to minimize the spin and orbital phase shift ac- 
quired by an incident Cooper pair—. Indeed, since the 
orbital phase change for the given heterostructure is at 
any angle of incidence of the Cooper pairs, the overall 
phase shift is minimized when the magnetization forms 
an angle tt/2 with the (i- vector— ~—. The weak modifi- 
cation of the pairing amplitude is related to the orbital 
symmetry of the triplet superconductor that is perpen- 
dicular to the interface. Very close to the interface in the 
TSC side there is also a small contribution coming from 
the inverse proximity effect which tends to favor the su- 
perconducting pairing amplitude when the ferromagnetic 
magnetization is in the (xy)-planc. 

On the other hand in the FM side the pairing ampli- 
tude exhibits a spatial damped oscillatory behavior at 
any angle 9 ^ tt/2 and a pure monotonous dependence 
when the magnetization is perpendicular to the d- vector. 
The oscillatory behaviour of the pairing amplitude comes 
from the spin pair-breaking effect of the ferromagnet. In- 
deed, the Cooper pairs spin is perpendicular to the d- 
vector, and therefore the only non pair-breaking config- 
uration for the superconducting components within the 
ferromagnet happens when 9 = tt/2. 

Moreover, the oscillation length does not depend 
on 9 because the split of the energy for the up and 
down electrons, that controls the modulation in the 
proximity, is always proportional to the amplitude of the 
ferromagnetic exchange h and not on the pair-breaking 
projection. It is worth pointing out that differently 
from the SSC junction, the oscillation of the pairing 
amplitude leads to a negative phase shift only for a value 
of 9 in the range ~ [0,7r/4]. Still, we notice that in this 
case there are no spin singlet induced components at 
the interface differently to what is observed in a planar 
FM-TSC hctcrostructures where the interface explicitly 
breaks the inversion symmetry in the spatial directions 
that enter into the spin-triplet order parameter—. 
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FIG. 2. (color online) Spatial variation of the real part of 
the p^-wave pairing amplitude F Vlc scaled to the TSC bulk 
value for various angles of the magnetization 8 with respect 
to the d-vector as defined in Fig. [TJ For this case we assume 
h = hi nt = 0.16 and 8 = p. The FM-TSC interface (i z = 0) 
is denoted by the dashed line. The inset shows the evolution 
of the pairing amplitude F Pm as a function of the angle 8 at a 
given layer i z — 2. 



FIG. 3. (color online) variation of the Gibbs free energy (solid 
black line) and the order-parameter derived energy E op (blue 
dotted line) as functions of the angle 8 at a given value of the 
ferromagnetic exchange h = 0.32. Both energies are scaled 
with respect to their minimum amplitude (AG = G — G mm 



and A E op = E op 



77 mill \ 
J^op I 



and are renormalized to the E 



value of a single superconducting layer in the homogenous 
state. 



We arc now interested in the energetically most fa- 
vorable magnetization orientation in the ferromagnetic 
layers. We have analyzed the evolution of the Gibbs free 
energy G and of the energy associated with the TSC order 
parameter E op for different values of the relative orien- 
tation of the ferromagnetic moment with respect to the 
d- vector and as a function of the strength of the ferromag- 
netic exchange. It is worth stressing that E op is related 
to the condensation energy through the amplitude of the 
density of states at the Fermi level for the single particle 
spectrum. Thus, if there occurs an electronic reconstruc- 
tion at the TSC layers close to the interface, due to the 
split Fermi surfaces in the fcrromagnct, this effect can 
modify the condensation energy and as a consequence it 
can in turn affect the total energy balance. 

As one can notice in Fig. |3]thc Gibbs free energy has a 
non monotonous profile with two degenerate minima at 
9 = ±7r/2 and a maximum for 9 = 0. This means that, 
due to the coupling with the TSC, the energetically most 
favorable magnetic configuration for the ferromagnet is 
achieved when the magnetization lies in the plane of the 
Cooper pairs spin. Otherwise, the maximal loss in the 
energy occurs when the magnetization is parallel to the 
d- vector. To understand the driving force in the stability 
of the magnetic configuration it is instructive to analyze 
the behavior of the order-parameter derived energy -E^p- 
The results indicate that the maximum energy gain due 
to the variation of the superconducting order parameter 
in the TSC region occurs for a magnetization that lies 
perpendicular to the d- vector which also corresponds to 
a lack of spin-pair breaking in the proximity effect. Since 



the total energy has the same minima of the E op term 
one can deduce that for this strength of the exchange field 
the changes of the spectra in the Gibbs free energy due 
to the proximity effect and to a possible electronic recon- 
struction at the interface cooperate with the gain in the 
condensation energy to stabilize the observed magnetic 
configuration. 



IV. SPIN-ACTIVE PROCESSES AT FM-TSC 
INTERFACE 

In this section we analyze the role of spin-dependent 
processes at the FM-TSC interface. Concerning the 
modification of the pairing amplitude, the presence of a 
spin-scattering interface does not modify substantially 
the trend observed for the uniform case with 9 = p 
where the TSC order parameter gets maximized when 
the magnetization lies in the (xy)-plane. This is mainly 
because the overall phase shift experienced by the 
Cooper pairs is again minimized when the magnetization 
is perpendicular to the d-vector. Nevertheless, other 
processes are induced due to the spin-active character of 
the interface such as triplet-singlet conversion and spin- 
filtering. These mechanisms contribute in determining 
the amplitude of the triplet pairing nearby the interface 
in the superconducting region of the heterostructure 
and the intensity of the proximity effect. The general 
outcome is that the pairing amplitude in the TSC 
region increases when 9 and p tend to be collinear 
and antiparallel while in the FM region it is larger for 
the parallel interface-bulk magnetic configuration if 
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compared to the antiparallel one. 



A. Optimal magnetic configuration 

The analysis of the energies variations reveals a com- 
plex scenario because the presence of spin-dependent 
processes at the interface tends to modify the energy 
contributions of the proximity effect to the total balance 
as well as induces a spin-dependent intcrfacial electronic 
reconstruction. In the following we consider two possible 
physical situations: i) the magnetization in the bulk 
layers and at the interface of the FM is isotropic and we 
look for the energetically most favorable magnetic state, 
ii) the magnetic moment of the ferromagnetic layers in 
the bulk is fixed in amplitude and orientation due to 
an intrinsic source of magnetic anisotropy or is pinned 
externally along a specific direction (for example via 
an exchange-bias coupling) and we look for the optimal 
orientation of the magnetization at the interface which 
minimizes the Gibbs free energy. For both cases the 
amplitude of the spin moments is chosen by fixing the 
value of the ferromagnetic exchange field h and h- ln t in 
the bulk FM layers and at the interface, respectively. 
The behavior of the Gibbs free energy G and the order- 
parameter derived energy E op with respect to the angles 
p and 9 are reported in Figs. 0] and 03 respectively, for 
two representative values of the ferromagnetic exchanges 
h = h int = 0.32 and 0.64. 

When considering the physical case i) where all the 
magnetic orientations are assumed to be equivalent 
within the FM, the presence of the spin-triplet super- 
conductor can lead to a breaking of the rotational sym- 
metry for both the magnetic state in the bulk and at the 
interface. In this circumstance, since 9 and p have to 
be taken as independent parameters, it is the absolute 
minimum of the Gibbs free energy to indicate the most 
favorable magnetic configuration for the FM-TSC system 
in the presence of spin-dependent processes at the inter- 
face. By inspection of Fig. 0] (a) one can observe that for 
h = hint = 0.32 the ground state corresponds to a config- 
uration with parallel interface-bulk magnetic moments, 
i.e. 9 = 7r/2 and p = tt/2, with the magnetization lying 
in the (xy)-plane perpendicular to the d-vector. How- 
ever, for h = h int = 0.64 (Fig. [5] (a)) the Gibbs free en- 
ergy is minimized by an antiparallel configuration within 
the (xy)-plane i.e., 9 = tt/2 and p = —tt/2. The behav- 
ior of the order-parameter derived energy shows that the 
maximum gain in E op always occurs for an antiparallel 
interface-bulk configuration. 

Moving to the case ii) where the magnetization in 
the FM bulk is kept fixed while it can rotate at the 
interface, it is possible to observe a different scenario 
for the optimal magnetic configuration of the FM-TSC 
heterostructure. Let us start from the exchange values 
h = hi n t = 0.32 (Fig. 0]). In general, the magnetization 
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FIG. 4. (color online) a) variation of the Gibbs free energy, 
scaled with respect to its minimum amplitude G mm (AG = 
G — G mln ) and renormalized to the energy E op associated with 
a single superconducting layer in the homogenous case, as a 
function of the angle p at various angle and for a given 
value of the ferromagnetic exchange h — h- m t = 0.32. b) 
Variation of E op (AE op = E op — E^ n ) energy associated 
with the superconducting order parameter in the TSC region 
with the same renormalization and parameters as in a). 



at the interface turns out to be non collinear to the one 
in the FM bulk if we exclude the in-plane configuration 
at 9 = 7r/2. Indeed, for a bulk magnetization paral- 
lel to the d-vector (i.e. 9 = 0) the optimal orientation 
at the interface is at an angle p « 3ir/8. The angle p 
corresponding to the relative minimum of the Gibbs free 
energy moves from about 37r/8 to 7r/2 when 9 varies from 
to 7r/2. Interestingly, the minimum of the order param- 
eter energy E op , corresponding to the maximal variation 
of the gap amplitude, occurs always at p = —tt/2 for any 
value of 9 if we exclude the 9 = configuration where 
±7r/2 are degenerate in energy due to symmetry reasons. 
Such behavior anticipates the scenario that the selection 
of a given orientation for the interface spin moment is 
not driven by an energy gain due to the superconduct- 
ing order parameter but rather it is a consequence of a 
mechanism of energy gain that originates in the FM re- 
gion due to the proximity effect. On the other hand, 
for h = h^t = 0.64 (Fig. O the minimum of the Gibbs 
free energy is always at about p = —tt/2 for any given 
orientation 9 of the bulk magnetization. 
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FIG. 6. (color online) phase diagram for the energetically 
most favorable configuration in the plane (hi nt ,h). The fer- 
romagnetic exchange fields are expressed in units of the ex- 
change amplitude Hum — 3.15 corresponding to the transition 
value to an half-metal FM, where namely the density of the 
spin minority carriers goes to zero. The sketch represents 
the most favorable configurations with parallel (white area) 
and anti-parallel (blue area) interface-bulk magnetic moments 
alignment, respectively. 



FIG. 5. (color online) a) variation of the Gibbs free energy, 
scaled with respect to its minimum amplitude G mln (AG = 
G — G mln ) and renormalized to the energy E op associated 
with a single superconducting layer in the homogenous case, 
as a function of the angle p at various angles 9 and for a 
given value of the ferromagnetic exchange h = hint = 0.64. 
b) Variation of E op (AE op = E op — E%p m ) energy associated 
with the superconducting order parameter in the TSC region 
with the same renormalization and parameters as in a). 



B. Phase Diagram 

The previous results indicate that the ground state is 
not unique but depends on the strength of the exchange 
field. Therefore it is interesting to calculate a phase dia- 
gram in terms of h and h- lIL t by evaluating the minimum 
of the Gibbs free energy in the [8, p] landscape. The 
results in Fig. [5] show that the parallel and antiparallel 
interface-bulk alignment of the magnetization lying on 
the (xy)-plane are the two most stable configurations. 
The first one is favored for a weak and strong ferromag- 
net while the second is preferred at intermediate coupling 
values. A first order transition occurs at the boundary 
between the two regions of stability. It is worth pointing 
out that there is an asymmetry in the role played by the 
interface and the bulk exchange field which can be em- 
phasized by scanning the phase diagram along lines with 
given h lnt and h. Indeed at small fixed values of h- ln t 
(< 0.2) there is no reentrant behavior as a function of h 
while one does find the double transition for the antipar- 
allel configuration in the case of small fixed values of h 



by tuning the amplitude of /i; nt . Such asymmetry under- 
lines the role of the interfacial spin-dependent scattering 
in driving the reentrant behavior from the parallel to the 
antiparallel configuration in the regime where the mag- 
netization at the interface is large with a small density 
of spin-minority electrons. 

To understand which mechanisms can lead to such a 
behavior, it is interesting to study the contribution of 
the superconducting order parameter energy E op as pre- 
viously shown. The Fig. [7] shows the energy difference 
between both configurations for hi nt = h. We can ob- 
serve that the gain in the condensation energy tends to 
favor the antiparallel pattern over the parallel one for 
any value of the interface-bulk ferromagnetic exchange. 
Hence, as a result of the comparison between the total 
Gibbs free energy and the order parameter energy E op , 
one can deduce that the competition between the an- 
tiparallel and the parallel configuration originates from a 
modification of the energy spectrum due to the proximity 
effect within the FM. 



C. Density of States 

To get a deeper understanding of the modification 
on the energy spectrum induced by the interface-bulk 
misalignment we evaluate the density of states (DoS) 
p(i t ,w) at each layer within the FM and TSC region 
of the heterostructure. We also analyze the integrated 
quantity I(i z ) = f_ D u)p(i z ,ui)duj, representing the 
contribution to the Gibbs energy coming from the 
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FIG. 7. (color online) evolution of the difference of the Gibbs 
free energy (black solid line) associated with the antiparallel 
and the parallel configurations (see sketch in Fig. [6]) as a func- 
tion of the exchange field strength h assuming that h = hi nt . 
The blue dashed line indicates the difference between E op of 
the antiparallel configuration and that one for the parallel 
magnetic state. 



FIG. 8. (color online) difference of the integrated density 
of states for the antiparallel and the parallel configuration 
as a function of the layer position and of the exchange field 
strength. We consider an equal amplitude of the exchange 
field at the interface and in the bulk, i.e. h = hint- 



electronic spectrum, with —D being the lowest occupied 
energy level and the effective Fermi energy for the 
system after the Bogoliubov rotation. 

We report in Fig. |8] the spatial evolution of the dif- 
ference between the values of I(i z ) corresponding to the 
antiparallel (i.e. p = ±7r/2 and 8 = z f7f/2) and the par- 
allel (i.e. p = 6 = ir/2) configurations as a function of 
the ferromagnetic exchange, assuming that h = /lint- 

We notice that I(i z ) is a negative quantity and thus 
a positive difference indicates that the electronic spec- 
trum contribution tends to favor the parallel configura- 
tion, while the antiparallel one is favored when the dif- 
ference is negative. The analysis of the DoS reveals a 
significant dependence of the contribution I(i z ) on the 
magnetic configuration and on the strength of the ferro- 
magnetic exchange. The main change in the DoS occurs 
at the FM layers which are close to the interface. This is 
a consequence of the interplay between the misalignment 
of the interface-bulk magnetizations, the scattering pro- 
cesses at the interface, and the leaking of the supercon- 
ducting pairs via the proximity effect. Since the interface 
electronic reconstruction tends to modify the density of 
states profile, it is expected to affect also the I(i z ) be- 
havior, thus playing a role in setting the most favorable 
magnetic configuration. 

The behavior of I(i z ) is strongly layer dependent. In- 
deed, as shown in Fig. [5J the difference between the am- 
plitude I(i z ) associated with the parallel and the antipar- 
allel configurations is vanishing away from the interface 
both in the FM and in TSC regions of the heterostruc- 
turc. Otherwise, it exhibits a nonmonotonous behavior 
for the closest FM layers to the interface when the fer- 
romagnetic exchange is varied from the weak coupling to 



the half-metallic regime. 

In the TSC region the electronic reconstruction always 
favors the antiparallel state, as also expected from the 
evaluation of the order parameter energy E op . Thus, no 
significant renormalization of the DoS close to the Fermi 
energy happens that interferes with the role of the order 
parameter in the energy balance. 

In the FM side the energy contribution due to the 
electronic spectrum is highly nontrivial. For small values 
of the ferromagnetic exchange there is an appreciable 
difference between the energy related to the electronic 
spectrum of the two magnetic configurations. Such 
difference becomes larger and larger in favor of the 
parallel configuration when h increases from 0.32 to 
0.51. Then the difference decreases until the antiparallel 
configuration starts to dominate. The results of the DoS 
confirm a scenario of electronic reconstruction close to 
the interface playing a relevant role in the configuration 
transition. Indeed, the nonmonotonous behaviour of 
the difference between the amplitude I(i z ) associated 
with the parallel and antiparallel configurations follows 
the same trend of the Gibbs free energy shown in Fig. 
[7] while the condensation energy gain contributes to 
cooperate for the stability of the antiparallel state. 

It is worth pointing out that the modifications of the 
DoS are substantial both in the energy range of the su- 
perconducting gap and also outside it. This is partially 
due to the superconducting correlations and in part to 
the change of the DoS due to the processes of scatter- 
ing at the interface at energies which are larger than the 
superconducting gap in the presence of exchange split 
Fermi surfaces. Interestingly, in the regime of intermedi- 
ate ferromagnetic exchange field, i.e. for h in the range 
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FIG. 9. (color online) Spatial variation of the real part of 
FJJ and pair correlations associated with the ft (solid 
lines) and 4-4- (dashed line) spin-triplet configurations for a 
representative case of p x orbital symmetry as a function of 
the magnetization direction 6 and at a given exchange field 
h — 0.16. The amplitude of the correlator is scaled to the 
value of the bulk one Re[F Px ] in the t4- spin configuration. 
The FM-TSC interface is denoted by the dashed line. 



of 0.32 to 0.57, the density of states of the ferromagnetic 
layers close to the interface exhibits new midgap elec- 
tronic states in the antiparallel configuration. Although 
such states may represent a distinctive feature of the FM- 
TSC heterostructure, they do not contribute significantly 
in determining the stability of the ferromagnetic profile. 

V. RE-ORIENTATION OF THE SPIN-TRIPLET 
PAIR CORRELATIONS 

In this section we analyze the induced spin-triplet pair 
correlations, as defined in Eq. |H with p x and p y or- 
bital symmetry but having parallel tt an d antiparallel 
4-4, spin configuration as compared to the state with zero 
spin projection associated with the bulk order parame- 
ter. These pair correlations are identically zero in the 
bulk, since the superconductor has no net spin polariza- 
tion. However, it is possible to induce equal spin pairing 
correlations close to the interface, since a magnetization 
noncollinear to the d- vector can lead to spin-flip processes 
and then to a consequent spin-mixing. It is instructive to 
consider the case of an uniform ferromagnet i.e., p = 9, 
with a magnetization forming an angle 9 with respect to 
the d- vector. 

In Fig. |H]we show the spatial profile of the F°°(i z ) pair 
correlations for both cases a = (t, 4-) for different orien- 
tations of the magnetization. As one can notice, there 
are no induced spin triplet pair correlations for a mag- 
netization parallel to the d- vector (9 = 0). By changing 
the magnetization orientation from the direction perpen- 
dicular to the (cn/)-plane (9 = 0) to in-plane (9 = n/2), 



the amplitude of the equal spin triplet correlations ex- 
hibits a nonmonotonous behavior. It increases until a 
critical angle 9 max ~ 7r/4 and then decreases when the 
magnetization moves toward the (xy)-plane. The evolu- 
tion of the pair correlations in the FM side is decaying 
and oscillating as a function of the layer position with 
the same period of those associated with the bulk order 
parameter (see for comparison the Fig. [5]). Since there 
is always a split between the Fermi energies of electrons 
with spin up and down at any orientation of the mag- 
netization, we do expect a spatial dependence with an 
oscillating behavior. Moreover, it is interesting to see in 
the TSC side that a pair correlation is also induced by 
inverse proximity effect. Therefore, in the case of the 
existence of an equal-spin pairing potential V" in such 
a system, there would be a rotation of the d- vector that 
can have significant consequences on the electronic trans- 
port in a Joscphson junction configuration^. Still, the 
modulation of the pair correlator does not lead to a phase 
change as a function of the spatial coordinate. 

A change in the amplitude of the ferromagnetic ex- 
change field does not modify the general trend obtained 
for the case of h = 0.16. In particular by increasing the 
exchange one gets a reduction of the equal spin triplet 
correlations. This is because the overall amplitude is set 
by the magnitude of the correlator associated with the 
superconducting order parameter. Concerning the spin- 
active processes, the interface-bulk misalignment tends 
to modify the amplitude of the equal spin pairing com- 
ponents close to the interface. 

It is worth pointing out that, since the maximum am- 
plitude of the equal-spin triplet correlations corresponds 
to an angle 9 max ~ 7r/4, they do not have a relevant role 
in the Gibbs free energy related to the proximity effect. 
On the other hand, in the case of an induced pair poten- 
tial in the TSC region in the equal-spin channel due to 
a modification of the pairing glue at the interface, there 
would be a significant modification of the energy balance 
related to the condensation energy. 



VI. CONCLUSIONS 

We have studied the physical behavior of a 
ferromagnet- spin-triplet superconductor heterostructure 
with a layered geometry assuming an orbital chiral sym- 
metry for the superconducting order parameter and spin 
dependent scattering at the interface. By means of a 
self-consistent approach on a three-dimensional lattice we 
have determined the spatial profile of the pairing ampli- 
tude and magnetization in terms of the strength of the 
ferromagnetic exchange and the spin-active processes at 
the interface. We find that the pairing amplitude in the 
TSC region close to the interface increases if the mag- 
netization moves from being parallel to perpendicular to 
the <i-vector. In the presence of a spin-active interface 
a configuration with in-plane antialigncd magnetic mo- 
ments for the bulk-interface layers maximizes the value 
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of the pairing amplitude. The penetration of the Cooper 
pairs in the FM leads to induced pairing amplitudes with 
zero-spin projection that we show have a damped oscil- 
lating behavior if the magnetization is out-of-plane, while 
it is monotonously decaying when the magnetization is 
perpendicular to the d- vector, i.e. parallel to the spin 
direction of the Cooper pairs. As a distinctive feature, 
the oscillating behavior brings a phase change only if the 
angle between the magnetization and the d-vector is in 
the ~ [0, 7r/4] window. The spin-active scattering mainly 
modifies the behavior close to the interface by influencing 
the intensity and the phase change of the pairing ampli- 
tude. It is worth pointing out that equal spin pairing 
correlations are also induced both at the interface in the 
TSC and with a long range penetration length in the 
FM if the magnetization is noncollinear to the <i-vector. 
We have shown that their amplitude is maximized at an 
angle 9 ~ 7r/4 for a uniform ferromagnet and it has a 
non trivial behavior in the presence of a bulk-interface 
magnetic misalignment. We have pointed out that such 
correlations can play a major role in the case of the exis- 
tence of a pair potential in the same spin channel at the 
interface of the spin-triplet superconductor. This effect 
can lead to a re-orientation of the (i-vector approaching 
the ferromagnet and drastically influence the anisotropy 
of the ferromagnet as well as the charge and spin trans- 
port of the hybrid structure^. 

The knowledge of the key physical quantities of the 
FM-TSC hybrid provides the basis for a detailed analy- 
sis of the optimal magnetic configuration as due to the 
coupling with a spin-triplet superconductor. The study 
of the Gibbs free energy of the system reveals that the 
rotational spin invariance of the ferromagnet is generally 
broken due to the coupling with the superconductor and 
in turn the FM develops an anisotropy which favors an 
in-plane orientation of the magnetization. In the case 
of a uniform ferromagnet, we have demonstrated that 
the driving mechanism for the stabilization of the bulk- 
interface relative orientation of magnetization is based 
on the gain of the superconducting condensation energy 
rather than on the proximity effect. The switching on of 
an active ferromagnetic layer at the interface leads to a 
different physical scenario. There are two optimal mag- 
netic states which minimize the Gibbs free energy of the 
system: they have a parallel or antiparallel bulk-interface 
alignment with the magnetization oriented in the plane. 
By exploring the full space of parameters as a function 
of the relative exchange field strength at the interface 
and in the bulk we have demonstrated that the system 
can switch between the two ground-state configurations. 



The antiparallel configuration is generally favored for a 
regime of weak and strong ferromagnet otherwise the 
parallel state is stabilized. In the latter the electronic 
reconstruction in the FM domain close to the interface 
competes with the energy contribution due to the mod- 
ification of the superconducting order parameter which 
would always tend to favor an antiparallel interface-bulk 
magnetic configuration. We emphasize that such elec- 
tronic reconstruction has important contributions also 
in an energy window outside the gap. Furthermore, we 
have demonstrated that the competition between the var- 
ious mechanisms involved in the interplay between ferro- 
magnetism and superconductivity in the FM-TSC hybrid 
leads to a noncollinear magnetic profile in the case of a 
strongly anisotropic bulk ferromagnet. 

The results obtained in this paper represent a solid ba- 
sis to investigate a FM-TSC hybrid in a conventional epi- 
taxially grown geometry. We point out that they can be 
of direct interest for the case of the Si~2Ru04 spin-triplet 
superconductor interfaced with other ferromagnetic ox- 
ides having matching lattice parameters, i.e. SrRu03 or 
doped LaMn03, but can also be indirectly used to extract 
information of the symmetry of the order parameter of 
unknown superconductors through the investigation of 
the magnetic profile in the ferromagnet. Indeed, the oc- 
currence of magnetization anisotropy within the FM can 
be experimentally probed measuring the magnetic mo- 
ment profile as a function of the layer position via nu- 
clear magnetic resonance, muon spin rotation or polar- 
ized neutron reflectometry. The latter has been success- 
fully applied, for instance, to observe a surprisingly large 
superconductivity-induced modulation of the ferromag- 
netic magnetization profile in superlattices made of high 
temperature superconductors interfaced with ferromag- 
netic manganitessi. Finally, since a spin-active interface 
plays a relevant role in determining the character of the 
spin-triplet pair correlations induced in the FM and in 
tuning the free energy profile, it is expected to be an im- 
portant ingredient for possible spin-valve devices based 
on spin-triplet superconductors as well as for controlling 
the dynamical evolution of the magnetization in the fer- 
romagnet. 
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